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We analyse the role of partonic transverse motion in unpolarized Semi-Inclusive Deep Inelastic 
Scattering (SIDIS) processes. Imposing appropriate kinematical conditions, we find some constraints 
which fix an upper limit to the range of allowed fcx values. We show that, applying these addi- 
tional requirements on the partonic kinematics, we obtain different results with respect to the usual 
phenomenological approach based on the Gaussian smearing with analytical integration over an 
unlimited range of k± values. These variations are particularly interesting for some observables, like 
the {cos (j)h) azimuthal modulation of the unpolarized SIDIS cross section or the average transverse 
momentum of the final, detected hadron. 

PACS numbers: 13.88. -|-e, 13.60.-r, 13.85.Ni 

INTRODUCTION 

Inclusive and Semi-Inclusive Deep Inelastic Scattering (DIS and SIDIS) processes are important tools to 
understand the structure of nucleons and nuclei. Spin asymmetries in polarized SIDIS are directly related to 
Transverse Momentum Dependent (TMD) parton distribution and fragmentation functions, and are the subject 
of intense theoretical and experimental studies. The usual, coUinear parton distribution functions depend on the 
fraction of hadron momentum carried by the scattering parton and on the virtuality of the probe, Q^. TMDs 
additionally depend on the intrinsic transverse momentum of the parton, fcj^, opening invaluable opportunities 
to unravel the three-dimensional partonic picture of the nucleon in momentum space. 

At leading twist, the spin structure of a spin- 1/2 hadron can be described by eight TMDs [T]-0]- TMDs 
represents particular physical aspects of spin-orbit correlations at the partonic level. The dependence of the 
SIDIS cross section on the azimuthal angle, (f)h^ of the electro-produced hadron with respect to the lepton 
scattering plane and on the nucleon polarization azimuthal angle, (jjs, allows a term by term separation of the 
different azimuthal contributions to the measured unpolarized and polarized cross sections and spin asymmetries. 

The unpolarized SIDIS cross-section can be used not only to study the unpolarized TMD distribution function 
fq/p{x, k±) and the unpolarized TMD fragmentation function Dfi/q(z,p±), that encode the intrinsic dynamics of 
unpolarized partons, but also the Boer-Mulders distribution and the Collins fragmentation functions, which carry 
information about the dynamics of transversely polarized partons inside hadrons and give rise, for instance, to 
a cos 2(f>h modulation of the unpolarized cross section. The existence of partonic intrinsic transverse momenta 
is also unequivocally signaled, in the unpolarized SIDIS cross section, by a cos (j^h modulation, which is a 
subleading twist effect, suppressed by one power of Q. This contribution to the unpolarized cross section 
consists of a purely kinematical term, the Cahn effect [U [S] , proportional to the convolution of unpolarized 
distribution and fragmentation functions, and a term proportional to the convolution of the Boer-Mulders and 
the Collins functions, together with other twist-3 contributions, as pointed out in Ref. [2]. 

Polarized SIDIS experiments, in addition, allow us to explore the other TMDs which describe the dynamics 
of polarized and unpolarized partons inside polarized nucleons, such as the helicity, transversity and Sivers 
functions. These will not be considered in this paper. 

TMD factorization theorems [B] are proven for SIDIS at momentum scales Q ^ Pt ^ ^qcd\ here large 
Q values are needed to allow a perturbative treatment of the underlying partonic subprocess, while the small 
scale Pt ensures that the observables are sensitive to the intrinsic parton motion. Moreover, the smallness of 
intrinsic quark momenta is also explicitly required, as this approach is based on a series expansion in terms 
of the ratio fcj^/Q and higher orders are usually neglected. However, the present knowledge of TMDs and 
azimuthal asymmetries is mainly based on the available experimental data from the HERMES and COMPASS 
Collaborations, that operate at relatively low values: typically, the cut > 1 GeV^ is assumed in order to 
ensure SIDIS kinematics. 

Very often, in phenomenological analysis, the transverse momentum distribution of the TMDs is assumed to 
be a Gaussian. Although very simple, this approximation leads to a successful description of many sets of data 
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(see for instance Refs. [71-fT7]V However, as the quality and the amount of data has significantly improved lately, 
a more detailed picture of the partonic description is now becoming more and more necessary. In particular, 
the latest COMPASS analysis [T8j of the SIDIS unpolarized cross section suggests that the average transverse 
momentum of the detected hadrons could sensitively depend on the parton momentum fractions xb and 
and even on flavor, although similar conclusions cannot be presently inferred from HERMES [1^1 dO] and 
JLAB data [H], as we will discuss in more details in Section [HI A In Ref. [S], a simple Gaussian model was 
used to successfully describe the Cahn cos (ph azimuthal modulation measured by the EMC Collaboration [22] 
and allowed the extraction of the unpolarized TMD's Gaussian widths. Later, it has been realized that the 
predictions of Ref. |9] largely overestimated the more recent data from HERMES US] and COMPASS [H] 
which have meanwhile become available. On the contrary, a large Cahn effect now seems to be required in 
order to describe the cos2(/>;i asymmetries [25], at least at the COMPASS kinematics. Furthermore, a global 
phenomenological analysis of data led by Schweitzer, Teckentrup and Metz [15j suggests that the average 
transverse momenta (A;^), used as a free parameter in the Gaussian, depends on y/s, the center of mass energy 
of the target and the incoming lepton. 

This complex scenario thus requires an accurate critical analysis of the work done so far, to try and shed 
some light on these apparent controversies. The detailed study of the COMPASS Collaboration on the Pt 
distribution of their SIDIS measurements '18], pointing to a considerable deviation from the expected behavior, 
prompted us to a reanalysis of the approximations used in our fc_L integration. In this paper, we examine the 
kinematical regions where order 0{k\/Q'^) and higher can be safely neglected and we realize that not all the 
kinematical domains of the existing and planned experiments satisfy the basic criteria k±^ ~ Pt — f^QCD Q- 
Inspired by the parton model, and still adopting a Gaussian model for our TMDs, we bound the integration 
range of transverse momenta and we observe, in some kinematical regions, remarkable deviations from the 
predictions obtained from the common TMD approach, based on the Gaussian parametrization integrated over 
the full k±^ range, [0,cx)]. We show that some kinematical ranges, typically low x^ or equivalently low 
regions, are not safely controlled by the present phenomenological model, while bounds on fcj_ can prevent 
uncontrolled large k^/Q contributions. This leads, for instance, to a better description of some observables like 
the (cos (f)h) asymmetry and introduces some interesting effects in the (P^) behaviors. However, simple parton 
model considerations are not sufficient to fully describe the present, in some sense puzzling, data: soft gluon 
emissions, higher twist contributions and QCD evolution of the TMDs can play a role. 

Higher order contributions to cross sections and asymmetries are difficult to estimate; here, we will consider 
only those coming from purely kinematical corrections to twist-2 TMDs and show that they are not negligible 
in most of the present experimental setups. This makes the extraction of twist-3 TMDs from existing data 
troublesome. 

Future experimental data from an Electron Ion Collider [211 HZ], where the range would be easily adjustable 
in order to estimate higher order contributions, will definitely help to disentangle them from leading twist 
contributions. 



I. SIDIS KINEMATICS 



Let us give a brief review of (.p — I'h + X Semi-Inclusive Deep Inelastic Scattering kinematics. We study 
the SIDIS process in the 7* — p cm. frame, where 7* denotes the virtual photon. Following the Trento 
conventions |28| . we take the virtual photon momentum q along the -\-z direction and the proton momentum 
P in the opposite direction, as shown in Fig. [T] The leptonic momenta define a plane that coincides with our 
x-z plane. The detected hadron has momentum Ph, its transverse component is denoted by Pt and (j)h is its 
azimuthal angle: together with z they identify the hadron production plane, Fig.[T] We adopt the usual SIDIS 
variables neglecting the lepton, the proton and the final hadron masses: 

s = {P + lf Q' = -q' (P + g)2 = I^2^^^^Q2 

Xb 

P-q 

"^ = 2p:^"W^W '^y-e^^s -^ = (P-P^)/(P-i)^ (1) 

where £ is the incoming lepton momentum. The proton and the virtual photon momenta can be written in the 
7* — p cm. frame, as functions of the invariants W and Q in this way: 

'^^U^-^'^'^'^+f) ^ = ^0(1,0,0,-1) Po-l(w+^)- (2) 
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FIG. 1: Kinematical configuration and conventions for SIDIS processes in the 7* — p cm. frame. The initial and 
final lepton momenta define the x-z plane, while the detected hadron momentum and the z axis determine the hadron 
production plane, at an angle (p^. 



In the parton model the virtual photon scatters off a on-shell quark. The quark momentum k can be written, 
in the 7* — p cm. as: 

^Po + -!^,k^,-xPo + S^) (3) 



4a;Pn 4xPn 



(4) 



where x — k~/P~ is the quark light-cone momentum fraction (see Appendix [A| for more details) and k± is 
the quark intrinsic transverse momentum, see Fig. [2] The final emitted quark has momentum k' — q -\- k. Its 
on-shell condition 

= 2g • fc - Q2 ^ (5) 

implies [7] 

(6) 




As indicated in Fig. [2] = Ph — {Ph ■ k ) k is the transverse momentum of the hadron h with respect to 
the direction k' of the fragmenting quark, while z — P^ /k'~^ is the light-cone fraction of the quark momentum 
carried by the resulting hadron. The fragmenting variables z and pj_ can be expressed in terms of the usual 
observed hadronic variables Ft and Zh, as shown in Eqs. (26) and (28) of Ref. [7]. 




FIG. 2: Kinematical configuration and conventions for the SIDIS partonic sub-process: k and fe' are the initial and 
final quark momenta, and fcj_ is their transverse component. Ph is the final, detected hadron momentum, with a 
component, transverse with respect to the fragmenting quark direction k' . 
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All the kinematical relations given above are exact expressions at all orders in a {kj_/Q) expansion. Neglecting 
terms of order 0{k'j_/Q'^), they considerably simplify and we find: 

x~xb, z~ Zh Pt - Zhkj_ . (7) 

Transverse Momentum Dependent distributions depend on the kinematical variables x and k± defined above. 
Let's consider, for instance, the unpolarized distribution function fq/p{x,k±), which gives the number density 
of unpolarized quarks inside an unpolarized proton; this function is usually normalized in such a way that 

d?kx_ fq/p{x, k±) = fq/p{x) , (8) 

where fq/p{x) is the usual, coUinear parton distribution function at some given scale Q^. The same logic holds 
for the unpolarized TMD fragmentation function, Dh/q{z,p±). 

Very often, in phenomenological analysis, a Gaussian dependence of the TMDs is assumed, adopting the 
following parametrizations: 

fq/p{x, k^) = fq/p{x) ^ e-^-i/<'=i> (9) 

and 

Dh/q{z,Pi_) - Dh/q{z) , (10) 

where fq/p{x) and D^/q^z) can be taken from the available fits of world data, while (fc^) and {p^) are free 
parameters to be extracted from dedicated fits. 

Note that the functions defined in Eqs. (|9]) and (10) obey the normalization condition of Eq. (|8|, so that 



27r /•CO 

dip / dk±_ kj_ fq/p{x, k±_) = fq/p{x) . (11) 
Jo 

The Gaussian parametrization in principle allows any value of kj_ from zero to infinity. However, the integrand is 
weighted by the Gaussian k± distribution, so that k± values larger than the Gaussian width, (fc^), are strongly 
suppressed. Typical values of {k'j_) are of a few hundreds MeV. Thus if is large with respect to (fc^), the 
Gaussian represents an effective model that prevents large k±/Q contributions to the cross section. However in 
many low energy SIDIS experiments, like HERMES and COMPASS, (Q^) ~ 2 GeV^ and the experimental cut 
> I GeV^ is used: thus at low x we have ~ 1 GeV^. Therefore, in these particular cases, the Gaussian 
smearing is not sufficient to cut away large k±/Q contributions. Notice that these considerations are still valid 
for any phenomenological parametrization that does not impose any cut to the range of allowed kj_ values. 
Hence the necessity to explore whether it is possible to find a physical picture that allows us to put some further 
constraints on the partonic intrinsic motion. 

II. PHYSICAL PARTONIC CUTS 

In order to find some constraints on the partonic intrinsic motion, we adopt the simple picture of nucleons 
provided by the parton model. Although "non-physical" in many aspects, the parton model can be seen as good 
approximation and a toy-model to understand some physical QCD features. In particular, it gives kinematical 
limits on the transverse momentum size, which can be obtained by requiring the energy of the parton to be 
less than the energy of the parent hadron and by preventing the parton to move backward with respect to the 
parent hadron direction (k^ < 0). The energy bound implies: 

xPo + ^ <Po kl< 4x(l - x)P^ 

^ ki< "i|'"\ Q^- (12) 
xb[^ - xb) 



Inserting Eq. ([6| in Eq. (12 1 and solving one finds: 

fci < (2 - xj(l - xJQ2 ^ < X, < 1 . (13) 
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phase space as determined by the bounds of Eqs. (13 1 and (15 I. 



FIG. 3: fci/Q^ 

both bounds, is represented by the shaded area below the solid line, corresponding to Eq. 
corresponding to Eq. ( |15[ ). Notice that present data from HERMES and COMPASS experiments span the region 
xb ^ 0.3, where only the momentum bound of Eq. ( 15 I plays a role. 



The allowed region, which fulfills 
( 13 1 and the dashed line, 



Requiring the parton to move in the forward direction with respect to the parent hadron gives: 

(P • fc) > ^ kl< p2 



xsi'^ - Xb) 



Using Eq. ([6| and solving we find 



(l-2xJ2 



(14) 



(15) 



Notice that these are exact relations, which hold at all orders in {k±/Q). These constraints are obtained 
assuming that partons are on-shell (the parton ofF-shellness should be very small, at least in the region where 

we will apply the parton model). 

The ratio k^/Q^, as constrained by Eqs. (13) and (15), is shown in Fig. [s] as a functions of x^: from this plot 
it is immediately evident that although in principle Eq. ( 15 ) (represent ed b y the dashed line) gives a stringent 
limit on fc^/Q^ in the region 



< 0.5, it intercepts the bound of Eq. ( |13| (solid line) in x^ ~ 0.3, where the 
latter becomes most relevant. Notice also that present data from HERMES and COMPASS experiments span 
the region xb ^ 0.3, where only the momentum bound of Eq. ( 15 ) plays a role. 

Once the maximum value of k± is bounded by the conditions of Eqs. (13) and (15), to comply with the 
normalization condition of Eq. (|8|) , we set the right normalization coefficient 



fq/p{x,kj_) fq/p{x) 



1 



-K/ikl) 



1 _ p-ikT'^V/ikl) TT{kl) 



(16) 



where (fc™'''^)^ denotes the maximum value of fc^ for each given values of and as required by Eqs. ( 13 ),( 15 ) 
so that 



fq/pix) 







dip I k±dk± fq/p{x,k±) . 
Jo 



(17) 



This normalization will allow us to reproduce correctly all inclusive cross-sections, such as the DIS cross- 
section, for which collinear results are known. Note that if Q — >■ oo then fc™"^ — >■ oo and one obtains the usual 
normalization of Eq. ([s]). 
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III. UNPOLARISED SIDIS CROSS SECTION 



According to Refs. [TH3] the unpolarized differential cross section for the SIDIS process, i + p - 
be written as 



£'hX can 



dXg dy dzh d'^Px 



^ { ^-^^^^-^ Fuu + (2 -y)y^ cos <j>,,F^°^^- + {l-y) cos2^^ F^°f^- } , (18) 



where the F "structure functions" , which involve the relevant convolutions of distribution and fragmentation 
functions over the intrinsic transverse momenta, are defined within a TMD factorization scheme, at order 
0(fci/Q), as 



Fuu — 



y d'^k_ifq/p{x,kj_)Dh/q{z,p_L), 



(19) 



pCOS (ph 



2E 



2 / ^2,^^ ^ 



e~ I d k 



{Pt ■ k±_) fg/p{x,kj_) Dh/giz,pj_ 
, Pt- zuk^{PT-k^) 



2pA 



Afgyp{x,kj_) A^D,^/g^(z,pj_) 



(20) 



/^^fc. 



Pt (Pt ■ fc±) - 2zhk^{PT ■ k^f + z^k^ 
2p± 



^fqyp{x, k±) A^ Dh/qf{z,Pl_) 



(21) 



Although three different azimuthal modulations are simultaneously at work in the total unpolarized SIDIS cross 



section of Eq. (18), to extract single effects one can introduce the following azimuthal moments 



(cos (ph) 



fdc^i.da^P^^'''^ 
/ dcj)h da^P 



COS( 
i'hX 



{cos2(j)h) 



I 



^'^^ cos20,j 



(22) 



which isolate the contributions of F'^^ '^'^ 



and F, 



cos 2</)h 



UIL. 



from the traditional, collinear term Fi 



uu- 



As we showed in 

(|19|, (|20|) and (|21|) coincide with those defined in Ref. considering 



Ref. P^, the F structure functions in Eqs. 
only leading twist TMDs. In the phenomenological study we will perform in what follows, we will not consider 
higher twist dynamical contributions [2 to the structure functions, therefore we shall keep in mind that strong 
deviations of measured quantities from the predictions obtained in this simple model will signal the presence of 
higher twis t co ntrib utio ns. 

In Eqs. (19) and (20), fq{x,k±) is the unpolarized TMD distribution function and Dl^{z,p±) is the unpo- 



larised TMU fragmentation function discussed in Section |Tj while A/^t/p(a;, A:_l) is the Boer-Mulders distribu- 
tion function, related to the number density of transversely polarized quarks inside an unpolarized proton, and 
Di^/qr{z,p±) is the Collins fragmentation function which, in turn, is related to the number density of trans- 
versely polarized quarks fragmenting into a spinless hadron. Other common notations used for the Boer-Mulders 
and Collins functions are: 



^fqyp{x,kj_) 



M 



hi {x, fcj 



A'^D 



h/q 



iz,p±) 



zMh 



Ht{z,pi_) 



(23) 
(24) 



where M and Mh are the masses of the initial proton and of the final hadron, respectively. 



Summarizing, Eq. ( 19 1 corresponds to the usual, collinear f^^p 
cross section; in the cos0/i azimuthal modulation of Eq. (20) the first term is the Calm effect [HIS], signaling 



D/j/q contribution to the unpolarized SIDIS 



in the cos ( 

the existence of transverse momenta in the partonic scattering even when considering only unpolarized partons, 
while the second term corresponds to the Boer-Mulders ® Collins contribution of transversely polarized quarks 
to the unpolarized cross section; finally, Eq. (21) gives the cos20/i azimuthal modulation of the unpolarized 
SIDIS cross section. An additional contribution to the cos (ph modulation can be generated by "dynamical" 
higher twist contributions from twist-3 functions [2]. Instead, the cos2</)/i azimuthal moment in Eq. (21) is 



not suppressed by 1/Q and does not receive any dynamical or kinematical twist-3 contribution [5], but can be 
affected by twist-4 contributions, like the cos2(l)h Cahn effect discussed in Section p^II C[ 
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Therefore from this aziniuthal moment one could extract the Boer-Mulders and Collins TMD's, in all those 
kinematical ranges in which higher orders in the k±/Q expansion can safely be neglected. As a matter of 
facts, as we will see in Section [nj higher orders of {k±/Q)"', n > 2, could potentially be important when 
analysing HERMES, COMPASS and JLab experimental data, and should not be neglected. Therefore twist-4 
contributions could be important in these kinematics. One of these contributions, the "twist-4" Calm effect, 
Eq. (351, has been analyzed phenomcnologically in Ref [IS] , see Section III C 

For the Boer-Mulders and Collins functions, we assume parametrizations similar to those in Eqs. ([9| and 
(10), with an extra multiplicative factor fcj^ or p±, respectively, to give them the appropriate behavior in the 
small fcj^ or p± region [21 [H] : 



^fqt/p{x,kj_) = Af t/p{x) V2e 



(25) 



N-. 



pj__ e 



-p1/{p1)c 



(26) 



with 



(kl) 



(pI) 



Ml 



(27) 



The x-dependent function A/^t/p(a;) and the z-dependent function A Z)^/^t(z) in Eqs. (25) and (26) are not 
known, and should be determined phenomcnologically by fitting the available data on azimuthal asymmetries 
and moments; the k±^ and pj_ dependent terms and their normalization are chosen in such a way that positivity 
bounds [2^ are fulfilled automatically. 

At this stage, the usual procedure would be to perform an analytical k±^ integration of Eqs. ( 19|pT| ) over the 
range [0, oo], using the parametrizations in Eqs. ([9]), (10), (25) and (26) with the appropriate normalization of 
Eq. (|8|, and to re-express all the F structure functions in terms of the Gaussian parameters: 



Fuu — 



pcos 2<j}h 



pCOS 0ft 



E 

9 



-Pt/{Pt)c 



9 



2 -g- E "^9 fq/pMDh/qizh) 



Zh{kl) 



-L/ BM 



(pI 



2 \2 



(28) 
(29) 



A/,t/p(a;J A^D^,/^,{z^) e-PT/{P^)nM 



(fci)L(pi ) 

{kl){pl 



,2 \2 , 



\ 11 BM 



^l{kl)sM{P^~{P^)s.)+{pl)c{PT) 



Tl BM 



(30) 



where 



{PT)BM={pl)c+^l{kl), 



(31) 
(32) 



We stress that the analytical integration which leads to Eqs. (28)-(30) is performed over the full range of fcj_ 
values: 



cPki 



dip / dk± k± 



(33) 



and that the expressions (|3l| and (32) which relate the (Px) to {p]_), {k^} and Zh are due to such choice. In 



other words, Eqs. ( |31[ ) and ( |32[ ) are a direct consequence not only of the assumption of Gaussian k± and p± 
distributions of the TMDs, but also of the choice of the k± integration range. Many phenomenological analysis 
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on the TMDs are based on the Gaussian assumption and overlook any issue regarding the Umits of integration 
over k±. However, great attention to this should be payed when analysing data from JLAB, HERMES and 
COMPASS experiments, where the average is not so large (~ 2 GeV^). Therefore in some kinematical 
ranges it could happen that the {k±/Q) values accessed are not small. One immediately visible signal of this is 
the Cahn effect in both azimuthal moments, at twist-3 in {cos (j>h) and at twist-4 in {cos2(j)h), which are directly 
proportional to k±/Q and k^/Q^ respectively, and are found to be (phcnomcnologically) large. 



A. Impact of the partonic cuts on the Fuu term of the unpolarized cross section and on 



In this section we will show the impact of the fcj^ cuts presented in Eqs. ( 13 ) and ( 15 ) on the calculation 



of the SIDIS unpolarized cross section and on the average transverse momenta of the final detected hadron, 
(P^). Figure [4] shows the first term of the SIDIS unpolarized cross section, proportional to Fjju, calculated 
for HERMES and COMPASS kinematics (the detailed experimental cuts are reported in Appendix |b]) for 7r+ 
production, integrated over all variables but one, x^, Zh and Pt- 

We use the unpolarized integrated PDF's given in Ref. 30 and the unpolarized fragmentation functions of 
Ref. |31| . For the Gaussian widths of the unpolarized distribution and fragmention TMDs we use the values 
extracted in Ref. [7j: (kj) = 0.25 GeV^ and (p^) = 0.20 GeV^. The solid (red) line, denoted in the legend 
as "num+cuts" , corresponds to the unpolarized cross section calculated according to Eqs. (18), (19|-(21) with 
a numerical kj_ integration over the range obtained implementing the fcj^-cuts of Eqs. (13 1 and (15). The 
dashed (blue) line, indicated as "analyf, is the unpolarized cross section calculated according to Eqs. (28)- 
(30), resulting from a k± analytical integration over the range [0, 00]. These plots clearly show that, as far as 
the and distributions are concerned, there is no difference between the two calculations. Instead, a slight 
modification can be observed in the Pt distribution (see the upper and lower right panels). 

The average hadronic transverse momentum (P^) of the final, detected hadron h is defined as: 



Jd^PTPTdcr 



(34) 



HERMES Proton-7i* 



num+cuts ■ 
analyt ■ 





0.2 0.4 0.6 0.8 




0.2 0.4 0.6 0.8 1 



COMPASS Proton-ji* 






0.2 0.4 0.6 



FIG. 4: Unpolarized cross section for tt^ production at HERMES (upper panel) and COMPASS experiment on a proton 
target (lower panel), as a function of (left plot), Zh (central p lot) and Pt (right plot). The (red) solid line corresponds 
to the unpolarized cross section calculated according to Eq. ( 18 1 with a numerical fcx integration implementing the 
fcx-cuts of Eqs. ( 13 1 and ( 15 1. The dashed (blue) line is the unpolarized cross section calculated according to Eqs. ( 28 1- 
(|30| resulting from an analytical kj 



integration from zero to infinity. We do not show the analogous cross section 
corresponding to the COMPASS experiment on deuteron target as the effects of the k± cuts are very similar to those 
for COMPASS on proton target. 
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HERMES Proton-Ti'^ 




0.25 
0.2 
0.15 



analyt - 
<pi)+2^(ki) - 



0.2 0.4 0.6 



COMPASS Proton-7i+ 



0.28 
0.26 
0.24 
0.22 
0.2 
0.18 
0.16 
0.001 



num+cuts - 
analyt - 



analyt - 
<pi>+2^(k^) - 



FIG. 5: (P|), defined in Eq. (34 1, as a function of (left plot) and of zf^ (right plot), for tt"*" production at HERMES 
(upper panel) and COMPASS (lo wer pan el). The solid (red) line corresponds to (Py) calculated with a numerical 
integration implementing Eqs. ( 13 1 and ( 15 1, while the dashed (blue) line is (Pt) calculated with an analytical integration. 
In both cases we have applied the experimental cuts on Pt reported in Appendix [b] . Finally, the dash-dotted (black) 
line corresponds to the Gaussian {Pt)g ■ 



Notice that if the integral in Eq. (34 1 is performed over the range [0, oo], then (P^) coincides with the Gaussian 
width of the unpolarized Pt distribution of Eq. (28): (P^) = {Pt)g- The experimental Pt range, however, 
usually span a finite region between some P™*" and PJp"^^] ther efor e, in a ny e xperimental analysis, one inevitably 
has {Pt) {Pt)g, even without considering the cuts in Eqs. (13) and ( |15| ). Consequently, the relation (P^) ~ 
(Pj.) + ^h(j^±) holds only approximatively. 

Figure shows the average hadronic transverse momentum (P^) as a function 
at HERMES and COMPASS, respectively. The sohd (red) lines correspond to 



of 



and of 



for 



( 18 ) with a numerical kj_ integrations and im 
(Pt) calculated according to Eqs 



Eq. 

correspond to 
taken into account the appropriate experimental 



calculated according to 
elementing Eqs. (|13[) and (|15[). Instead, the dashed (blue) lines 



P|) 



(|28|)-(30) with an analytical integration. In both cases we have 
cuts on Pt reported in Appendix [B| Clearly, at low x, there is 



as a function of zf^. 



a substantial deviation from the analytical calculation, which also affects the value of 
As far as the dependence is concerned, first of all, one can see that there is a large deviation from the naive 
formula, Eq. (31), corresponding to the dash-dotted (black) lines, for both calculations. Secondly, although the 
z^-dependence is not linear any more, it seems to be approaching an almost constant behavior (as opposed to 
what COMPASS latest analysis (TS] seems to indicate). Notice that, independently from Eqs. (13) and (15), if 
we naively assume (P|.) = (p^) -I- Zh{k\) in any analysis of the data we would conclude that the effectivefA:^) 
is smaller than the value (/c^) = 0.25 GeV^ that we used as an input in the Gaussian. Also (p^) ^ 0.20 GeV^, 
as a consequence of the limits of integration on Pt- Fig. [6] shows how the integratio n ra nge in flue nces the value 
of {Pt) when we integrate analytically (left plot) or implementing the cuts in Eqs. (13) and (15) (right plot). 

Finally, different behaviors in could, in principle, imply different values of (P^)Tor 7r+ and tt^. We have 
explored this possibility in order to account for the slight discrepancy in the {Pt) corresponding to tt"*" and 7r~ 
as observed by COMPASS. Unfortunately the difference predicted by our model for 7r+ and 7r~ is extremely 
tiny, even at COMPASS kinematics, as can be appreciated in Fig. [Tj One should keep in mind that the tiny 
separation observed here occurs for purely kinematical effects, although larger differences between the 7r+ and 
7r~ average Pt could be generated by adopting different k_\_ distribution widths for different quark flavours. 

In Fig. [8]we show how a cut in can change the description of data. We can see that cutting at higher 
means cutting the lower Xg region, where the constraint of Eq. ( 15 ) strongly applies. As a consequence, if we 
apply high cuts the description of the data with or without fcj^-cuts is the same. This means that, provided 
Eqs. (13) and (15) are right, there is a region of or where we can safely assume that {k±/Q) corrections 
are small or negligible and where a phenomenological analysis is safe and unambiguous. 
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FIG. 6: {P|), at COMPASS kinematics as a function of zf^ for different ranges in Pt, calculated analytically (left plot) 
and numerically implementing Eqs. (131 and (151 (right plot). 
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FIG. 7; {Pt) as function of zl for 7r+ and vr" production at HERMES (left plot) and COMPASS (right plot) pion 
production on a proton target. 
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B. Impact of the partonic cuts on the azimuthal moment (cos^h) 
At this stage, we are ready to evaluate the effect of the physical partonic cuts on the (cost/)^) and (cos2^;i) 



azimuthal moments, Eq. (22), which represent the most delicate terms of the SIDIS unpolarized cross section. 



The {cos 4>h) modulation receives two contributions, both suppressed by one power of {k±/Q), see Eq. (20). 
The Cahn term, which is proportional to the convolution of the unpolarized distribution and fragmentation 
functions, was extensively studied in Ref. [7|- There, EMC measurements [22 on the cos (fih modulation and 
of the Pt distribution on the unpolarized SIDIS cross section were used to determine the Gaussian width of 
the k± distribution of the unpolarized distribution function fq/p{x, kj_). The second term is proportional to the 
convolution of the Boer-Mulders distribution function and the Collins fragmentation function and was neglected 
in Ref. [?]• Since then, new and higher statistics experimental data have become available [231 [53]: it is therefore 
timely and interesting to evaluate its net contribution to the {cos (j)h) azimuthal moment, using some reasonable 
estimate of the Boer-Mulders and Collins TMDs from the literature [131 [55] . 

Figure [9] shows how a large deviation from the analytical integration results is obtained by applying the k± 



bounds of Eqs. (13) and (15) when computing the Cahn effect contribution to {cos (ffh) corresponding to the 
HERMES and COMPASS kinematics, whereas for the EMC set up, see Fig. [TT] one can hardly notice any 
difference. The reason is quite obvious: while EMC used high > 5 GeV^ cuts, HERMES and COMPASS 
experiments typically have > 1 GeV^. 

In Fig. |9] our results, obtained with and without k± - cuts, are compared to the latest HERMES [23] and 
COMPASS [23] data. Although still showing a considerable deviation from the experimental data, our calcula- 
tion confirms that physical partonic cuts have a quite dramatic effect in the small x region, and should therefore 
be taken into account in any further analysis of these experimental data. 

To evaluate the influence of the partonic cuts on the contribution to (cos(/)h) originating from the Boer- 
Mulders(8) Collins term, we use the parametrization of Ref. [32] for the Collins function while for the Boer-Mulders 



function we apply the extraction of Ref [53]. It can be seen from Fig. 10 that the Boer-Mulders contribution 



is very tiny (it gives a correction of a few percents over the dominant Cahn contribution) and is not strongly 



affected by kinematical cuts of Eqs. (13) and (15 



The residual discrepancy between the model prediction and the measurements of the (cos 0^) azimuthal 
moment could indicate that higher twist contributions, from pure twist-3 functions, for example, might be non 
negligible in this modulation. More elaborated phenomenological studies including twist-3 TMDs would be 
necessary to confirm these observation. 



C. Impact of the partonic cuts on the azimuthal moment (cos20h) 



The (cos 2(f)) azimuthal modulation, at twist-2, consists of a term proportional to the Boer-Mulders(8) Collins, 



see Eq. (21 1. Although it is not affected by any twist-three corrections, in kinematical ranges where Q is 



not very large compared to the average kj_, twist-4 contributions cannot be neglected. In particular, a twist-4 
"Cahn-like" effect actually gives a large contribution to the (cos 20) azimuthal moment, as was explained in 



details in Ref. [5S]. This provides an additional term to the F^'^^'^'" structure function, Eq. (|21[), of the form 



2{Pt ■ k^f - 1 /,/p(x, k^) Dt,,q{z,p^) . 



(35) 



In Ref. [5S] , the presence of a twist-4 term proved to be crucial to understand the available experimental data 
from HERMES [331 133] and COMPASS [33 [3g results. Here, a detailed study inspired by the fact that the 
HERMES Pt spectrum can be reproduced by Monte Carlo calculations with (fc^) = 0.18 GeV^, pointed out 
that slightly different values of TMD widths might be required for different experiments. 



We revisit these calculations by applying the physical cuts on the partonic transverse momenta, Eqs. ( 13 1 and 



( 15 ), and find that the value of the twist-4 Cahn effect is very sensitive to the k± constraint, as it can be seen in 



the upper panels of Fig. [121 '^i^ii'^ the Boer-Mulders contribution does not exhibit such a strong dependence, as 
shown in the lower panels of Fig. 12 This is explained by the fact that, indeed, the accessed values of (fc^/Q^) 
are strongly suppressed by limiting the range of k_\_. The sum of the Boer-Mulders and Cahn-like contributions 
which reproduces the (cos 2<f>h) azimuthal modulation is presented in Fig. 13 One can see that the description 
of the available data is very good. In Ref. J5j a satisfactory description of the HERMES data was achieved by 
adopting a smaller value of {k\) = 0.18 GeV^ for the HERMES data, while keeping {k\) = 0.25 GeV^ for fitting 
the COMPASS data (see FIT II). Here a similar improvement is achieved by taking into account the physical 
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FIG. 9: Cahn contribution to {cos for tt"*" and 7r~ production at HERMES and COMPASS kinematics, as a function 
of a:^ (left plot), Zh (central plot) and Pt (right plot). The solid (red) line corresponds to (cos^h) calculated according to 
Eq. \2Q\ w ith a numerical fcx integration over the range [0, fc™"^]. The dashed (blue) line is {cos 4>h) calculated according 
to Eq. ( |30[) o btained by integrating over k±_ analytically. We do not show the Boer-Mulders contribution as it is negligible 
(see Fig.|10p. The fuU circles are prehminary experimental data from HERMES [53] and COMPASS [H] Collaborations. 



cuts on the partonic transverse momenta which, by cutting the range of allowed k±^ values, effectively reduces 
the average value of fcj^ decreasing the contribution generated by the Cahn effect. Nevertheless, a slight puzzle 
still remains: while HERMES data seem to demand a very small Cahn contribution, it can be seen from the 
analysis in Ref. [2S] and from Fig. 15 that COMPASS data seem to require a large Cahn contribution. Large 
Cahn contributions can only be generated by large k_\_ values, as we have seen. Since the COMPASS target is 
not a pure Deuterium target, but a Li^D target, possible nuclear effects can enhance the values of fcj^. Partonic 
transverse motion generated by the nuclear smearing mechanisms does not have to fulfill the bounds in Eq. ( 13 1 
and ( 15 ) and can be effectively simulated by a traditional Gaussian smearing, without any restriction. Future 
COMPASS data on pure hydrogen target will help our understanding. 
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FIG. 11: On the left and the central plots, respectively, the (Pt) and the (cos (f>h) azimuthal modu lation for tt"'" production 
at EMC. The solid (red) line corresponds to the calculation performed starting from Eq. (181 and then integrating it 
numerically, implementing Eqs. ( |13[ ) and ( |15[ ). The dashed (blue) line correspond to the analytical calculation. On 
the right plot the Cahn (red solid line) and the Boer-Mulder (blu dashed line) contributions to the (cos^h) azimuthal 
modulation. 
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IV. CONCLUSIONS 

In this paper we have studied SIDIS processes within a QCD parton model in the framework of TMD 
factorization; the dependence on the parton instrinsic transverse momentum is modeUed through a Gaussian 
parametrization. By requiring the parton energy to be smaller than the energy of its parent hadron and 
preventing the parton to move backwards relatively to its parent hadron, we were able to determine an upper 
bound, k™°-^ {xg , Q'^) , to the range of allowed values of k±. Under these assumptions, we then re-calculated 
the three terms of the unpolarized SIDIS cross section {Fuij, F^^'^'^ and F^*^^"^'*), and the detected hadron 
average transverse momentum, (P^). Notice that we made sure that the unpolarized TMD distribution function, 
fq/pi-'^B I respected the proper normalization condition by requiring that the integral over k± in the restricted 
range [0, ^"""^(2;^ , Q^)] would still give the usual, collinear fq/p{xg). 

Although the effects of our - cuts over the azimuthal-independent term Fujj are almost irrelevant (only a 
slight difference in the Fun dependence on Ft can be appreciated), we realized that the detected hadron average 
transverse momentum, (Py), and the azimuthal moments {cos(f)fi) and (cos20/j), are strongly sensitive to the 
constraints on the k± allowed values. In particular, by limiting the k± integration range, which effectively 
reduces the Gaussian width (fc^) suppressing the asymmetry at low a;^, (and consequently low Q^) values, 
we improve the description of the {cos cf>h) azimuthal moment data from HERMES [23] and COMPASS [2l] 
Collaborations, which were largely overestimated by the predictions obtained with an analytical k± integration 
over the unlimited k± range. Although the overall size of (cos0/i), a factor 2 (ore more) smaller than that 
obtained without k± - cuts, is in good agreement with the most recent experimental data, some discrepancies 
remain in the shape of our predictions, which can be interpreted as a signal of the existence of higher twist 
contributions, that have been neglected in our study. 

As far as the (cos 2(j)h) azimuthal moment is concerned, the situation remains slightly unresolved. Infact, while 
the effective reduction of the Gaussian width (fc^) helps in obtaining a satisfactory agreement with HERMES 
measurements [M;, some considerable inconsistencies remain in the description of the COMPASS data [37] 
which, instead, seem to suggest the presence of a much larger k^/Q^ Cahn contribution, that could only be 
achieved by increasing the average k±. This might suggest the presence of nuclear smearing effects in the Li^D 
target. COMPASS future data on a proton target will help to clarify this issue. 

Finally, we observed a significant deviation of the detected hadron average transverse momentum, (P^), from 
the theoretical value {P^)q = {p\) + zf^{k\), obtained by an analytical k± integration over the unlimited k± 
range. This is induced by two different mechanisms: on one side, the constraints we applied on the fcj^ range of 
integration and, on the other side, the (inevitable) Pt cuts operated in the experimental analysis. In general, 
we predict a much flatter behaviour than that of {P^)g, but yet quite far from that suggested by the COMPASS 
analysis of Ref. |18| . This is indeed an issue to be further studied in future, as only very preliminary data are 
presently available from COMPASS, HERMES and JLab. 

In this study, higher twist contributions were neglected, together with different mechanisms to generate the 
intrinsic transverse momenta, like soft gluon emission and TMD QCD evolution effects fSS] I3S]- Therefore 
more refined phenomenological descriptions are required to fix the details of such complex kinematics and 
dynamics. However, we have shown that some extra care should be taken when dealing with the present 
available experimental data, as they span a kinematical region in which k±/Q contributions can be large and 
are not safely under control, unless some limiting prescription over the allowed values of k± is applied. 

Future experiments, like the Electron Ion Collider (EIC) [2BJ[37], where the experimental cuts and the 
range would be easily adjustable in order to avoid unsafe kinematical regions, will definitely help us to gain 
a much clearer understanding of the three-dimensional structure of hadrons and, in particular, to disentangle 
higher order corrections from leading twist contributions. 
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Appendix A: Sudakov decomposition 

For the treatment of the SIDIS kinematics, we use the usual Sudakov decomposition for four vectors: 



where Hghtcone vectors and are 



• 7i_ = 1 



^ (1,0,0,1) , = -^(1,0,0,-1), 



(Al) 

(A2) 
(A3) 



Using hghtcone coordinates a 
as 



± 



l/-\/2(a° ± a^), and notations v^^ — [v ,v±\ we can rewrite these vectors 

(A4) 



<=[0-,l+,Oj = [l-,0+,Oj ,4 = (0,0,t;j.) 
The product of two four vectors is then 



V ■ W — V^W + V — v± ■ w± 



(A5) 



Note that vl^ ■ vt^ = s-nd n± ■ vt — 0. The momentum four-vectors corresponding to the proton, the 
virtual photon and the struck quark are: 



2P- "+ ' 

^71 _ -I 



2^ 



(A6) 



xP n1 



2xP 



where the Sudakov vectors , n_ are defined in Appendix [Aj a; = fc^ /P^ is the quark light-cone momentum 
fraction and fc_L is the quark intrinsic transverse momentum vector, see Fig. [l] Note that, according to the 
Trento conventions [28], P~ is the "large" component of the proton's momentum, i.e. along nt- 
The variables of Eqs. ( A6) then can be expressed as 



V2W ' 



P- 



1 

71 



w 



w 



(A7) 



and thus the proton and the virtual photon momenta can be written in the 7* — p cm. frame, as functions of 
the invariants W and Q, in this way: 



\/2W' V2' 



0, 



1 



w ^ 



w 



o+,o^ 



w 

P",0^,-P0) , 



where 



(A8) 
(A9) 

(AlO) 



In the Generalised Parton Model the virtual photon scatters off an on-shell quark, thus we will neglect the 
virtuality of the quark fc^ and write the quark momentum in the 7* — p cm. frame as: 



k"^ 

xP-, — ^,fej 
' 2a;P- ' 



i-2 1.2 



(All) 
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Appendix B: Kinematical cuts 



If not stated differently in the text, in our analysis we adopt the following kinematical cuts. For the HERMES 
experiment: 

g2 > 1 GeV^ , >10 GeV^ , 0.05 < Pt < 1.0 GeV (Bl) 

0.023 < a;^ < 0.27 , 0.2 < z,, < 0.75 , 0.3 < y < 0.85 , 

and for COMPASS: 

Q2 > 1 GeV^ , >25 GeV^ , 0.1 < Pt < 1.0 GeV (B2) 

0.003 <x^ < 0.13 , 0.2 < Zh< 0.85 , 0.2<y < 0.9 . (B3) 

Notice that hese kinematics correspond to the experimental cuts performed by these collaborations in their 
most recent analysis of the unpolarized azimuthal asymmetries [231 137) . 
For the EMC kinematics we used the following kinematical cuts [22]: 

g2 > 5 GeV^ , <90 GeV^ , 0.01 < Pt < 2.0 GeV (B4) 

Eh >5GeY, 0.1 <zh< 0.9, 0.2 < y < 0.8 . (B5) 

The Electron Ion Collider (EIC) is a proposed facility |26l [27] to provide further information on the proton 
structures. Thus we used the following hypothetical kinematical configuration 

Q2 > 1 or g2 > 10 GeV^ , > 2b GeV^ , 0.05 < Pt < 1 GeV (B6) 

0.2 < zh < 0.8 , 0.05 < 2/ < 0.8 . (B7) 



[1] P. J. Mulders and R. D. Tangerman, Nucl. Phys. B461, 197 (1996). 

[2] A. Bacchetta et al, JHEP 02, 093 (2007), hep-ph/0611265. 

[3] M. Anselmino et al, Phys. Rev. D83, 114019 (2011), arXiv:1101.1011. 

[4] R. N. Cahn, Phys. Lett. B78, 269 (1978). 

[5] R. N. Cahn, Phys. Rev. D40, 3107 (1989). 

[6] X. Ji, J.-P. Ma, and F. Yuan, Phys. Rev. D71, 034005 (2005), hep-ph/0404183. 

[7] M. Anselmino et al, Phys. Rev. D71, 074006 (2005), hep-ph/0501196. 

[8] M. Anselmino et al, Phys. Rev. D72, 094007 (2005), hep-ph/0507181. 

[9] M. Anselmino et al, Phys. Rev. D73, 014020 (2006), hep-ph/0509035. 
[10] M. Anselmino, M. Boglione, A. Prokudin, and C. Turk, Eur. Phys. J. A31, 373 (2007), hep-ph/0606286. 
[11] M. Anselmino et al, J. Phys. G36, 015007 (2009), arXiv:0710.1569. 
[12] M. Anselmino et al, Eur. Phys. J. A39, 89 (2009), arXiv:0805.2677. 
[13] M. Anselmino et al, Phys. Rev. D75, 054032 (2007), hep-ph/0701006. 
[14] M. Anselmino et al, Phys. Rev. D79, 054010 (2009), arXiv:0901.3078. 

[15] P. Schweitzer, T. Teckentrup, and A. Metz, Phys. Rev. D81, 094019 (2010), arXiv:1003.2190. 

[16] H. Avakian et al, CLAS Collaboration, (2010), arXiv: 1003.4549. 

[17] M. Boglione and S. Melis, (2011), arXiv:1103.2084. 

[18] J.-F. Rajotte, COMPASS Collaboration, (2010), arXiv: 1008.5125. 

[19] A. Airapetian et al, HERMES Collaboration, Phys. Lett. B684, 114 (2010), arXiv;0906.2478. 
[20] F. Giordano, (2008), DESY-THESIS-2008-030. 

[21] M. Osipenko et al, CLAS Collaboration, Phys. Rev. D80, 032004 (2009), arXiv:0809.1153. 
[22] J. Ashman et al, European Muon Collaboration, Z. Phys. C52, 361 (1991). 
[23] F. Giordano and R. Lamb, HERMES Collaboration, PoS DIS2010, 106 (2010). 
[24] G. Sbrizzai, COMPASS Collaboration, (2009), arXiv:0902.0578. 

[25] V. Barone, S. Melis, and A. Prokudin, Phys. Rev. D81, 114026 (2010), arXiv:0912.5194. 

[26] A. Deshpande, R. Milner, R. Venugopalan, and W. Vogelsang, Ann. Rev. Nucl. Part. Sci. 55, 165 (2005), hep- 
ph/0506148. 

[27] T. Horn, P. Nadel-Turonski, and C. Weiss, (2009), arXiv:0908.1999. 

[28] A. Bacchetta, U. DAlesio, M. Diehl, and C. A. Miller, Phys. Rev. D70, 117504 (2004), hep-ph/0410050. 
[29] A. Bacchetta, M. Boghone, A. Henneman, and P. Mulders, Phys.Rev.Lett. 85, 712 (2000), hep-ph/9912490. 
[30] M. Cluck, E. Reya, and A. Vogt, Eur. Phys. J. C5, 461 (1998), hep-ph/9806404. 



19 



[31] D. de Florian, R. Sassot, and M. Stratmanri, Phys. Rev. D75, 114010 (2007), hcp-ph/0703242. 

[32] M. Anselmino et ai, Nucl. Phys. Proc. Suppl. 191, 98 (2009), arXiv:0812.4366. 

[33] R. Lamb and F. Giordano, HERMES Collaboration, Nucl. Phys. A827, 225c (2009). 

[34] F. Giordano and R. Lamb, HERMES Collaboration, AIP Conf. Proc. 1149, 423 (2009), arXiv:0901.2438. 

[35] W. Kafer, COMPASS, (2008), arXiv:0808.0114. 

[36] A. Bressan, COMPASS, (2009), arXiv:0907.5511. 

[37] G. Sbrizzai, COMPASS Collaboration, (2010), arXiv: 1012.4910. 

[38] S. M. Aybat and T. C. Rogers, (2011), arXiv: 1101. 5057. 

[39] J. Collins, "Foundations of Perturbative QCD" , Cambridge University Press, Cambridge (2011). 



